Nonlinear shallow ocean wave soliton interactions on flat beaches 
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Ocean waves are complex and often turbulent. While most ocean wave interactions are essentially linear, 
sometimes two or more waves interact in a nonlinear way. For example, two or more waves can interact and 
yield waves that are much taller than the sum of the original wave heights. Most of these nonlinear interactions 
look like an X or a Y or two connected Ys; at other times, several lines appear on each side of the interaction 
region. It was thought that such nonlinear interactions are rare events: they are not. Here we report that such 
nonlinear interactions occur every day, close to low tide, on two flat beaches that are about 2,000 km apart. 
These interactions are closely related to the analytic, soliton solutions of a widely studied multi-dimensional 
nonlinear wave equation. On a much larger scale, tsunami waves can merge in similar ways. 
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The study of water waves has a long and storied history, 
with many important applications including naval architec- 
ture, oil exploration, and tsunami propagation. The mathe- 
matics of these waves is difficult because the underlying equa- 
tions are strongly nonlinear and have a free boundary where 
water meets air; there is no comprehensive theory. Here we 
report that X, Y, and more complex nonlinear interactions fre- 
quently occur on two widely separated fiat beaches and are 
not rare events, as was previously thought. In fact, these X- 
and Y-type interactions can be seen daily, shortly before and 
after low tide. These phenomena are closely related to the an- 
alytical solution of a multi-dimensional nonlinear wave equa- 
tion that has been studied extensively since 1970 HIE) and 
is a generalization of an equation studied by Korteweg and 
de Vries in 1895 [3|, which gave rise to the concept of soli- 
tons |4). From the universality of the underlying equation [5 1 
and the fundamental nature of these waves, it is expected that 
similar X- and Y-type structures will be seen in many different 
physical problems, including fluid dynamics, nonlinear optics, 
and plasma physics. 



BACKGROUND AND INTRODUCTION 

Water waves have been studied by mathematicians, physi- 
cists, and engineers for hundreds of years. While there are 
many types of water waves, here we will discuss solitary 
waves in shallow water; they are often called solitons and they 
have unique properties. Solitary waves in fluids [6 1 and oceans 
1 7 1 are a major and active research area. 

J. S. Russell, a naval architect, made the first recorded ob- 
servation of a solitary wave in the Union Canal, Edinburgh in 
1834: a stopping barge set off a solitary wave that went along 
the canal for one or two miles without changing its speed or its 
shape [8 1. He did experiments and found, among other things, 
that the wave's speed depends on its height; so he concluded 
that it must be a nonlinear effect. J. Boussinesq [9] in the 
1870s and D. Korteweg and his student G. de Vries in 1895 
derived approximate nonlinear equations for shallow water 
waves. They found both solitary and periodic nonlinear wave 
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solutions to these equations; they also found that the speed is 
proportional to its amplitude — bigger waves move faster. So 
Russell's observations were quantitatively confirmed. 

Between 1895 and 1960, solitary waves were mostly stud- 
ied by water wave scientists, mathematicians, and coastal en- 
gineers. In the 1960s, applied mathematicians developed ro- 
bust approximation techniques and found that the Korteweg- 
de Vries (KdV) equation appears universally when there is 
weak quadratic nonlinearity and weak dispersion Q. In 1965, 
Zabusky and Kruskal [4 | found that the solitary waves of the 
KdV equation have remarkable elastic interaction properties 
and termed them solitons. Gardner, Greene, Kruskal, and 
Miura [10] then developed a method for solving the KdV 
equation with rapidly decaying initial data; this method has 
been extended to many other nonlinear equations and is called 
the Inverse Scattering Transform (1ST) IfTTllTZll — such equa- 
tions are called integrable. 

In 1970, Kadomtsev and Petviashvili [1| (KP) extended 
the KdV equation to include transverse effects; this multi- 
dimensional equation, like the KdV equation, is integrable 
[2j. Our observations in this article are related to soliton solu- 
tions of the KP equation that do not decay at large distances; 
these interacting, multi -dimensional line soliton solutions can 
be found analytically [ 1 1 1. Before our observations, there was 
only one well-known photograph of an interacting line soliton 
in the ocean and it was thought that such interactions are rare 
events; it was taken in the 1970s in Oregon (Fig. 4.7b in ref 
ifTTln and is similar to Fig. [3] Since the KP equation has other 
X, Y, and more complex line soliton solutions, we sought and 
found ocean waves with similar behavior (Figs.[T]-[6]l. Surpris- 
ingly, these X, Y, and more complex types of line solitons ap- 
pear frequently in shallow water on two relatively flat beaches, 
some 2,000 km apart! These freely propagating, interacting 
line solitons are remarkably robust. While these interactions 
are not stationary and so only last a few seconds, a casual ob- 
server will be able to see them with the insights provided in 
this article. Interestingly, in laboratory experiments involving 
internal waves emanating from the interaction of cylindrical 
wave fronts, Maxworthy [131 Fig. 11] reported an X-type in- 
ternal wave interaction; Weidman et al. [ 14 1 later showed that 
the length of the stem in lf]~3l Fig. 1 1] follows a Hopf bifurca- 
tion when plotted against the intersection angle. 



2 




1 1 ■"»-* .• 




FIG. 1. A plot and a photograph of an X-type interaction, (a) A plot 
of an analytical line-soliton interaction solution of the KP equation 
at t = using {3| and (|4j. In this and the following plots, we picked 
the kj and Ps to be qualitatively similar to the photograph in part (b). 
Here, h = k 2 = 1/2, Pi = -P 2 = 2/3 so e A < 2 » 2.3. (b) Taken in 
Mexico on 3 1 December 20 1 1 ; notice the large amplitude of the short 
stem, (c) A 3d-plot of the solution in (a), which qualitatively agrees 
with (b); we only include one 3d-plot because the density plots show 
the interaction behavior clearly. 



OBSERVATIONS 

Single line, solitary water waves are familiar to every beach 
goer: they are localized in the direction of propagation and 
have a distinctive, hump-like wave profile. These waves break 
when they are sufficiently large compared to the depth and 
they often curve from transverse beach and bottom effects. 
We will focus on interacting line solitary waves that form X, 
Y, and more complex interactions. 

It was thought that X-type ocean wave interactions hap- 
pen infrequently. This is not the case: X- and Y-type ocean 
wave interactions occur daily, shortly before and after low tide 
on relatively flat beaches. M.J. A. observed these interactions 
near 20°41'22"N, 105°17'44"W in Nuevo Vallarta, Mexico 
from 2009 to 2012 between December and April. D.B. ob- 
served these interactions near 33°57'52"N, 118°27'35"W on 
Venice Beach, California in May 2012 — about 2,000 km 
away. Figs.[T]-[6]shows a few of the thousands of photographs 




FIG. 2. A plot and photographs of a Y-type interaction, (a) ki = 1 /2, 
k 2 = 1, P { = 3/4, P 2 = 1/4 so e Ail - = 0. (b) Taken in Mexico on 6 
January 2010. (c) Taken in California on 3 May 2012. 



that we took. The water depth where we saw these interac- 
tions was shallow, usually between 5 and 20 cm; the beaches 
are long and relatively flat; the interactions usually happen 
within 2 hours before and after low tide; the cross-waves pro- 
duced near a jetty appear to help induce these interactions. 
We found that these X- and Y-type interactions usually come 
in groups, which last a few minutes. We saw many X- and Y- 
type interactions each day that we made observations; the rel- 
ative frequencies of the interactions were different at the two 
beaches — M.J. A. saw X-type interactions like Fig. [T]more 
often than D.B. We also saw more complex interactions, such 
as three line solitons on one side of the interaction region and 
two line solitons to the other side, which we will call a 3-in-2- 
out interaction; these more complex interactions are much less 
frequent than X- and Y-type interactions. Our observations in- 
dicate that X- and Y-type interactions are remarkably robust: 
they typically persist through bottom-depth changes, pertur- 
bations from wind and spray, and sometimes even breaking! 

We observed three types of X interactions: an interaction 
with a short stem (Fig. [T}, an interaction with a long stem 
where the stem height is higher than the incoming line solitons 
(Figs. [5] and Q, and an interaction with a long stem where 
the stem height is lower than the tallest incoming line soliton 
(Fig. [5J. The amplitude of the short-stem X-type interaction 



FIG. 3. A plot and photographs of an X-type interaction with a 
longer stem, (a) k x = k 2 = 1/2, P a = -1/4 - 1CT 2 , P 2 = 3/4 so 
e An ss 51. (b) Taken in California on 2 May 2012 in shallower water 
than Fig-fTJ). (c) Taken in California on 4 May 2012. 



FIG. 4. A plot and photographs of an X-type interaction with a very 
long stem, (a) k t = k 2 = 1/2, P 2 = -P x + 10~ 10 = 1/2 so e A "- « 
5 x 10 9 . (b) Taken in Mexico on 28 December 2011 in shallower 
water than Fig.[5]5. (c) Taken in California on 3 May 2012. 



MATHEMATICAL DESCRIPTION 



The KP equation (TJ, 



can be quite large in deeper water. Interestingly, the length of 
the stem often increases as the depth decreases. Fig.|2]shows a 
typical Y-type interaction. A more complex interaction, with 
three 'incoming' and two 'outgoing' segments, is shown in 
Fig-© 



When one knows what to look for and when and where to 
look for them, X- and Y-type interactions are fairly easy to ob- 
serve. In addition to happening less frequently, more complex 
interactions are harder to see because they are highly non- 
stationary and have shorter interaction times than X- and Y- 
type interactions. Another difficulty is that most water waves 
break before X- or Y-type interactions form; so sustained ob- 
servation may be needed. Along with the photographs here, 
we have also taken many videos that show the development 
and general dynamics of these waves; the readers can watch 
some of these videos and see many more photographs at our 
websites fl5l . 
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is the two-space and one-time dimensional equation that gov- 
erns unidirectional, maximally-balanced, weakly-nonlinear 
shallow water waves with weak transverse variation. Here, 
sub-scripts denote partial derivatives, rj = rj(x, y, t) is the 
wave height above the constant mean height h, g is gravity, 
y=l — t/3, t = T lipgh 2 ) is a dimensionless surface tension 
coefficient, and p is density. When there is no y-dependence, 
the equation reduces to the KdV equation [3 |. The KP equa- 
tion was first derived in the context of plasma physics Q] and 
was later derived in water waves lfl6l . The sign of y is im- 
portant: there is 'large' surface tension when y < and this 
equation is called KPI; there is 'small' surface tension when 
y > and this equation is called KPII. We can rescale ([TJ into 
the non-dimensional form [5] 



(t^ + 6uu x + u xxx ) x + 3cruyy — 0, 



(2) 
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FIG. 5. A plot and photographs of an X-type interaction, where the 
stem has a lower rather than a higher amplitude, (a) fc[ = l,fo = 1/2, 
Pi = 1/2 - 1(T 7 , P 2 = so e A < 2 « 5 x 1(T 8 . (b) and (c) Taken in 
California on 3 May 2012. 



where u relates to the wave height 77 and cr — +1 corresponds 
to the sign of y. 

For large surface tension, KPI has a lump-type solution that 
decays in both x and y but has not yet been observed. Only 
recently has a large-surface-tension one-dimensional soliton 
been observed IfTTl : it satisfies the KdV equation and is a de- 
pression from the mean height. 

We will only discuss KPII here because surface tension is 
small for ocean waves. The KPII equation has solutions with 
a single-phase, which we will call line-solitons. We are inter- 
ested in the interactions of line solitons. These solutions can 
be found by so-called direct methods IfTTl : special /V-soliton 
solutions of the KP equation can be written in the form 1 1 8 1 



u — un — 2 



d 2 F A 



(3) 



dx 2 

where is a polynomial in terms of suitable exponentials 
This solution is convenient for finding the simplest such solu- 
tion: the first three are 



Fx = 1 + e ni , F 2 = 1 + ■ 



f 3 - l + \^ e 1* + V e 1i+1j+Aij +e i7i+ift+iB+Ai2+Ai3+A 2 3 ) 

1 </<3 1 <;< /<3 

(4a) 



FIG. 6. A plot and photographs of a 3-in-2-out interaction, where 
there are three line solitons on one side of the interaction region and 
two line solitons on the other side, (a) k { = 1, k 2 = 2, k 3 = 3, 
Px = -1/3, P 2 = -2/3, P 3 = -5/3. (b) and (c) Taken in California 
on 4 May 2012. 



where rjj = kj[x + Pjy - (kj + 3crPpt] + rjj\ kj, Pj, rjj' are 
constants, and 



(kj - kj) 2 - QjPj - Pj f 
(k . +kj) 2_ (r( p._p j) 2 



KJ. 



(4b) 



For KPII (where cr = 1), u\, Fx corresponds to the sim- 
plest one line soliton, which is essentially one-dimensional. 
The more interesting case of u%, F2 corresponds to the in- 
teraction of two line soliton waves. These interactions have 
distinct patterns: when e A ' 2 = 0(1), we get an X-type inter- 
action with a short stem (Fig. [TJ; when e A ' 2 » 1, we get an 
X-type interaction with a long stem where the stem height is 
higher than the incoming line solitons (Figs. |3]and|4|; when 
e A ' 2 <k 1, we get an X-type interaction with a long stem where 
the stem height is less than the height of the tallest incoming 
line soliton (Fig. [5}; and when e A ' 2 = 0, we get a Y-type inter- 
action (Fig. [2]). As mentioned earlier, the length of the stem 
appears to be correlated to the depth of the water. Short stems 
where e A ' 2 = 0(1) are usually found in much deeper water 
than long stem X- or Y-type wave interactions where e A ' 2 » 1 
ore A ' 2 <K 1. 

Recently, novel and exotic web-like structures for the KP 
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equation (A^-in-M-out) have been found using Wronskian 
methods [ 19 , 20 1 that go beyond the simplest 'building block' 
solutions of X- and Y-type line soliton solutions. Note also 
that an AMn-M-out solution (where M < N) can be found 
by starting with F# and taking kj and Pj such that e Au - N = 
■ ■■ = e AN -'- N = 0; Fig. [6] shows such a 3-in-2-out interaction. 
It was recently shown that these line interactions persist under 
the next order perturbations in the equations for water waves 
ETl : while the stem can be four times the height of the incom- 
ing line solitons in the KP equation, it is less than four times 
the height when higher-order terms are included. 



Il26l that significantly amplified the tsunami and its destructive 
power. This effect is remarkably similar to the initial forma- 
tion of an X- or Y-type wave: while it is initially a linear super- 
position effect, the interaction can be significantly modified 
or enhanced by nonlinearity after propagating to shore. More- 
over, for large distances (in the open ocean direction) an earth- 
quake induced tsunami will propagate approximately like the 
KP equation. So strong nonlinear effects from X- or Y-type in- 
teractions can have serious effects for land much further away; 
the destruction in Sri Lanka from the 2004 Sumatra-Andaman 
earthquake induced tsunami is an example of such a long dis- 
tance effect. 



X- AND Y-TYPE STRUCTURES AND TSUNAMI 
PROPAGATION 

Miles 11221 |231 first discovered that Y-type solutions could 
be associated with the KP equation; he also related it to 
"Mach-stem reflection", the phenomenon that occurs in gas 
dynamics. Interestingly, Wiegel E4ll reported that the 1946 
Aleutian earthquake induced tsunami caused a Mach-stem re- 
flection along the cliffs of the western edge of Hilo Bay in 
Hawaii. Yeh et al. l25ll revisited Mach-stem reflection in wa- 
ter waves with an inclined bottom, both analytically in the 
context of the KP equation and in a laboratory water wave 
tank. 

Recent observations of the 2011 Japanese Tohoku-Oki 
earthquake induced tsunami indicate that there was a 'merg- 
ing' phenomenon from a cylindrical-wave-type interaction 



CONCLUSION 

We reported that X- and Y-type shallow water wave interac- 
tions on a flat beach are frequent, not rare, events. Casual ob- 
servers can see these fundamental wave structures once they 
know what to look for. Extensive ocean observations reported 
here enhance and complement laboratory and analytical find- 
ings. We expect that similar interactions will be observed in 
many other fields — including fluid dynamics, nonlinear op- 
tics, and plasma physics — because the leading-order equa- 
tion here is also the leading-order equation for many other 
physical phenomena. 
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